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ABSTRACT 

Neptune Trojans and Plutinos are two sub-populations of Trans-Neptunian Objects located in the 1:1 and the 3:2 mean motion 
resonances with Neptune, respectively, and therefore protected form close encounters with the planet. However, the orbits of these 
two kinds of objects may cross very often, allowing a higher collisional rate between them than with other kinds of Trans-Neptunian 
Objects and a consequent and size distribution alteration of the two sub-populations. 

Observational colors and absolute magnitudes of Neptune Trojans and Plutinos show that: i) there are no intrinsically bright (large) 
Plutinos at small inclinations; ii) there is an apparent excess of blue and intrinsically faint (small) Plutinos; and iii) Neptune Trojans 
possess the same blue colors as Plutinos within the same (estimated) size range do. 

For the present sub-populations we analyze the most favorable conditions for close encounters / collisions to occur and address if there 
could be a link between those encounters and the sizes and/or colors of Plutinos and Neptune Trojans. We also perform a simultaneous 
numerical simulation of the outer Solar System over 1 Gyr for all these bodies in order to estimate their collisional rate. 
We conclude that orbital overlap between Neptune Trojans and Plutinos is favored for Plutinos with high libration amplitudes, high 
eccentricities and low inclinations. Additionally, with the assumption that the collisions can be disruptive creating smaller objects not 
necessarily with similar colors, the present high concentration of small Plutinos at low inclinations can thus be a consequence of a 
collisional interaction with Neptune Trojans and such hypothesis should be further analyzed. 

Key words. Kuiper Belt - Celestial mechanics - Methods: N-body simulations - Methods: data analysis 



1. Introduction 

Trans-Neptunian Objects (TNOs), also known as Kuiper Belt 
Objects (KBOs), are a population of small and primitive icy bod- 
ies orbiting (mostly) beyond Neptune. Their study is one of the 
most significant ways to obtain information on the early ages of 
the Solar System. 

Based on some distinct dynamical properties, the TNOs can 
be subdivided in several different sub-populations, often also 
called families or groups. Subdivide them as a function of their 
physical properties seems to be far more complex. TNOs have 
surface colors so diverse that can go from blue/neutral (i.e. solar- 
like) to extremely red. A possible explana tion for the w i de va- 
riety of colors was originally proposed by iLuu & Jewittl (l996) 
with the collisional resurfacing model. In this model, the com- 
petition between surface reddening, due to cosmic -ray bombard- 
ment, and a resurfacing with frozen material (assumed to be 
bluer) withdrawn from beneath its crust by impact collisions 
could be responsible for the observed wide range of surface col- 
ors. As the model seems to fail, at least in its simple form, more 
complex forms have been proposed. Na mely, collisional resur - 
facing combined with cometary activity (Delsa nti et al.l l2004). 
and collisional resurfaci ng with layered reddening and even re- 
bluing by cosmic -rays (Gil-Hutton 2003). An alternate idea pro- 
poses that surface colors are primordial (Tegler et al. 2003b, 
and references therein). More recently, Grundvl d2009l) shows 
that an object could lose its redness simply by ice sublimation 



(without resurfacing). Such could explain the lack of red col- 
ors among TNOs that evolved into Jupiter family comets due 
to their closeness to the Sun. Yet, as the aforementioned work 
acknowledges, the existence of the blue/neutral TNOs at helio- 
centric distances where ices do not sublimate suggests the ex- 
istence of other coloring mechanisms. Our understanding on the 
origin and eventual alteration of TNOs colors is still very limited 
and, up to the present, none of these approaches lead to a fully 
consistent explanation for t he color diversity (for a review see 
iDoressoundiram et"ai1l2008l) . 



Among the TNOs, two sub-populations caught our attention: 
the Neptune Trojans and the Plutinos. Neptune Trojans are the 
small bodies trapped in a 1 : 1 mean motion orbital resonance with 
Neptune, also with orbital eccentricities lower than 0. 1 . Roughly 
they co-orbit with Neptun e concentrated 60° ahead a nd 60° be- 
hind Neptune's position. Chiang & Lithwick (2005) proposed 
that Neptune Trojans formed by in-situ accretion from small- 
sized debris once Neptune's migration has stopped. On the other 
hand, other works indicate they can survive planetar y migration 
dNesvornv & Donesl l2002t iKortenkamp et al]|2004l) and, more 
recently, due to the disco very of a large population of high incli- 
nation Neptune Trojans. [Nesvornv & V okrouhlickv (2009) sus- 
tain they were captured during planetary migration Compared to 
TNOs Neptune Trojans seem to be quite small in size (diameters 
D < 100 km) and also slightly blue. 
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At the time of our analysis 6 Neptune Trojans were knowrfl 
all of them librating around the Lagrangian point L4. 

Plutinos are those trapped in a 3:2 mean motion orbital 
resonance with Neptune, with eccentricity values ranging be- 
tween 0.1 and 0.3. Unlike Trojans, the colors of Plutinos vary 
from blue/neutral to the very red, and their sizes range from 
a few tens of km to a few thousands (note that Pluto is a 
Plutino). Roughly, Plutinos possess semi-major axes within 39 < 
a < 40.5 AU. Through long-term dynamical evolution studies 
iLvkawka & Mukail d2007l) identified 98 Plutinos and that will be 
our reference. 

Being locked at the 3:2 resonance, Plutinos can periodically 
cross the orbit of Neptune without colliding with it. However, 
this protection from collisions is not possible for Neptune 
Trojans and Plutinos. A first look at the geometry of these two 
populations suggests that they might even collide frequently. 

The analysis of the collisio nal resurfacing mod el by 
IThebault & Doressoundiraml d2003l) and lThebaulU 120031) found 
that Plutinos were significantly more affected by collisions 
than other TNOs. Since, observationally, Plutinos did not ex- 
hibit bluer colors than all other TNOs taken as a whole, 
the aforementioned work strongly argued against the colli- 
sional resurfacing model, at least as major cause for the color 
diversity of TNOs. Neptun e Trojans were not included in 
IThebault & Doressoundiraml s simulations, though, since they 
were not known at the time. Notwithstanding the arguments 
against it, the collisional resurfacing that was analyzed was a 
very simple model. For instance, the possibility of collisional 
disruption of differentiated bodies producing fragments or ruble- 
piles with surface colors distinct from those of the parent bod- 
ies has not been taken into account in the collisional resurfac- 
ing models yet. Naturally, the possible outcomes of this scenario 
would be far more complex and our understanding of the physics 
of collisional processes between th e icy TNOs is still limited 
(see review bv lLeinhardt et al.l l2008). Nonetheless, both a colli- 
sional family of TNOs associated with (136108) Haumea and a 
collisiona l system of satelli t es associated with Pluto have been 
detected dBrown et al.ll2007tlStern et al.ll2006tlCanupll2005l) . 

Surveys show a power-law for the size distribution o f TNOs 
with slope change at D ~ 100 km dBernstein et alJl2004h . which 
is most consistent with TNOs being gravity dominated bodies 
with negligible material strength. Objects larger than ~ 100 km 
are difficult to disrupt, hence likely primordial, and objects 
smaller than this are expected to be shattered due to disru ptive 
col lisional evolu ti on dKenyon et alJl200l iPan & SariH"2005h . 

Ide Ella et alJ d2008h studied the collisional evolution of 
Plutinos considering only Plutino-Plutino collisions. They infer 
that any eventual slope change at D ~ 40 - 80 km on the power- 
law size distribution of Plutinos should be primordial and not 
of collisional origin. Yet, they find that collisional populations 
may form from the breakup of objects larger than 100 km, and 
if those populations form at low inclinations their fragments will 
likely stay in the resonance. 

On Sect. |2] we will analyze the observational data relative 
to Neptune Trojans and Plutinos. We will discuss two partic- 
ular observational properties: 1) there is an apparent "excess" 
of small blue Plutinos and these happen to be also in the same 
size range and to possess the same colors as the known Neptune 
Trojans, and 2) there is a concentration of small Plutinos, and 
absence of large ones, at low orbital inclinations. The eventual 
collisional interaction between Plutinos and Neptune Trojans 
has not been analyzed in previous works. Since their geometry 



points to the existence of mutual collisions, the existence of an 
eventual link between the colors and sizes of Neptune Trojans 
and Plutinos and their interactions/collisions motivated us to 
this work. On Sect. [3] we will study the dynamics of Neptune 
Trojans and Plutinos, separately, focusing on their orbits around 
the Lagrangian point L4. On Sect. [4] we will discuss the possi- 
bility of collisions between Neptune Trojans and Plutinos and 
the best conditions for this to happen. Finally, on Sect. [5] we will 
present our conclusions. 

2. Observational Results and Discussion 



In Table I A. 1 1 we summarize the orbital elements, B - R colors, 
and R-filter absolute magnitude s (Hr) for 4 Neptune Trojans 
(from Sheppard & Truiill ol2006l) and 41 Plutinos (see references 
on the table). In Fig. Q] we plot the orbital inclination vs orbital 
eccentricity of all these objects, together with their B-R colors, 
indexed on a color palette on the right hand side of the figure, 
and in which objects are plotted proportionally to their estimated 
diameter also indexed on the top of the figure. 
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Fig. 1: Orbital inclination vs eccentricity, estimated size, and color of 
Neptune Trojans and Plutinos for which those properties have been 
measured. 



The diameters of these objects (in kilometer s) are estimate d 
from the absolute magnitudes using the formula (Russell 1916): 



D, 



(km) 



2x 



2.24 X 10 16 x lo°-4(-2v.io-H fl ) 



PR 



(1) 



See http://cfa-www.harvard.edu/iau/lists/NeptuneTrojans.html 



where Hr is the R-filter abso lute magnitude, and the R-filter 
albedo is taken as p R = 0.09 dBrown & Truiilloll2004b . A size 
estimation exception is made for Pluto which is represented with 
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D — 2390 km. A vertical dashed line separates Neptune Trojans, 
located on the left side of the figure, from Plutinos, located on 
the right side. 

A first look at Fig.Q]shows that: 

(i) all Trojans are blufl small (D < 100 km) compared to the 
size distribution of Plutinos, and possess small eccentricities; 

(ii) there is an apparent concentration of small Plutinos at low 
inclination values and a concentration of large Plutinos at 
high inclinations; 

(iii) Plutinos have larger eccentricity values than Trojans, and 
their colors range from blue to red appearing randomly dis- 
tributed in inclination and eccentricity. 

(iv) all Plutinos within the same (estimated) size range as 
Neptune Trojans possess blue colors. 

It is important to note that the small sized Plutinos, which 
are also as blue as Neptune Trojans, are randomly scattered in 
eccentricity and inclination whereas the low inclined Plutinos 
are all relatively small but range from blue to red colors. 

Based on these properties we decide to explore two possible 
scenarios: 

(1) could the equally blue colors of the equally sized Plutinos 
and Neptune Trojans be the result of some collisional inter- 
action between both populations? 

(2) could the concentration of small Plutinos at low inclinations 
be the result of some collisional interaction between them 
and Neptune Trojans? 

Why these two scenarios? Fro m the simulati ons by 
iThebault & Doressoundiraml ( l2003l) and lThebaulil d2003l) . which 
analyzed the collision rates among all TNOs simultaneously 
(except for Neptune Trojans that were not known at the time), 
Plutinos received more collisions than any other sub-population. 
The number/energy of those collisions also seemed independent 
of their orbital parameters. That is, Plutinos seemed under the 
same collisional environment regardless of their high or low in- 
clination and/or eccentricity values. 

Therefore, since the possible effects of a significant number 
of Neptune Trojans were not considered in the aforementioned 
simulations in this work we explore if and how these Trojans 
could be the major cause of the non-homogeneity of surface 
properties observed among Plutinos and the (apparent) homo- 
geneity of surface properties observed for Neptune Trojans. 

Let us grasp scenario (1). The members of the collisional 
families associ ated with (136108) Haumea possess very similar 
bluish colors dBrown et alj 12007). It is reasonable to hypothe- 
size that if Neptune Trojans collided heavily with part of the 
Plutino population the collision outcomes could be similar in 
color. Considering that assumption, for this scenario to be possi- 
ble we need to find a similar collision rate between Trojans and 
Plutinos independent of the eccentricity and/or inclination values 
of Plutinos since the small (D < 100 km) and blue Plutinos are 
homogeneously scattered both in eccentricity and inclination. 

As to scenario (2), we are considering a less strict hypothesis 
in which the collision outcomes of collisions between Neptune 
Trojans and part of the Plutino population are not necessarily 
equal in color. The Pluto system has, presumably, a collisional 
origin, however, while Charon, Ni x, and Hydra possess similar 
colors Pluto is quite distinct (e.g. ISternll2009l) . Let us assume 



that the disruption of large layered TNOs could generate several 
objects with distinct colors either from fragments coming from 
different layers with different colors or from reaccretion of dif- 
ferent parts of a heterogenous cloud of debris. Hence, collisions 
would generate both blue and red, small and medium objects 
(D < 300 km). For this scenario to be possible we need to find 
a much higher collision rate between Trojans and the Plutinos 
with low inclination values than with the Plutinos at high incli- 
nations. Also, given that we do not see neither red nor mid/large- 
sized Neptune Trojans, whereas we see it among Plutinos in the 
presumable region of collisions, the primordial size and color 
distributions of Trojans and Plutinos should not have been simi- 
lar. 

Note that even though resonant objects are known to pos- 
sess l arge variations of their e ccentricity and/or inclination val- 
ues, iNesvornv & Roigl ([2000) showed the existence of two dy- 
namical populations of Plutinos, separating at i ~ 10°, which 
do not seem to easily mix. Consequently, any different colli- 
sional/surface evolution that both populations might have had 
in the past should still be detectable today. 

In order to test if the low inclined Plutinos are indeed 
small er than the other ones we use the Kolmogorov-Sm irnov 
test dPress etal.lll99H lKolmogorovlll933l: ISmirnov|[T939h . We 
successively compare the absolute magnitude Hr distribution of 
our 41 Plutinos (with measured B - R colors) having i < i c with 
those having i > i c , varying i c from 4° to 15° in increments of 
0.5°. Two significant solutions exist: a) Plutinos with i < 8° are 
intrinsically fainter, and assumed to be smaller, than those above 
8° with a (one-tailed) significance p = 0.0024; b) Plutinos with 
i < 13° are intrinsically fainter than those above 13° with a sig- 
nificance p = 0.0030. Even though we do not obtain an unique 
solution for the best inclination value that separates the larger 
(brighter) from the smaller (fainter) Plutinos, we can conclude 
that the low-inclined ones are in fact smaller. Figure [2] shows the 
cumulative distribution functions of absolute magnitudes for the 
higher inclined and lower inclined Plutinos in both solutions. 

We will proceed with the study of the dynamics of Plutinos 
and Neptune Trojans and investigate their possible collisions. 



3. Dynamical Evolution 

In previous section we discussed the possible relations between 
the colors of Plutinos and Neptune Trojans, and the eventual col- 
lisions between them. It is now our goal to test the possibility 
of collisions numerically. For that purpose we will simulate the 
outer Solar System evolution, where the TNOs are considered 
massless. This hypothesis is essential to speed up the integra- 
tions. The equation of motion for planets is given by 



i*i + Q(m s + m,) 



r, 



2> 



Yj - r, 



v r, 



'I, 



(2) 



2 For simplicity throughout this work we will call an object blue 
when B - R < 1.5 and red when B - R > 1.5. The B - R color of 
the Sun is 1.03. 



where r,- is the vector position of the planet, Q the gravitational 
constant, m s the mass of the Sun, m, the mass of the planet, and 
Np is the total number of planets. In our simulations we will only 
take into account the four giant planets. The effect of the inner 
Solar System in the dynamics of the Kuiper belt objects is only 
residual and by neglecting it we may use a larger step-size for 
numerical simulations and considerably decrease the length of 
the simulations. 
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Fig. 2: Cumulative distribution function of absolute magnitudes H R for 
Plutinos above and below 8° of inclination (top) and Plutinos above 
and below 13° (bottom). The incompatibility of distributions between 
the higher and lower inclined Plutinos is similar in both cases. There is 
an evident lack of intrinsically bright (large) objects among the lower 
inclined Plutinos. 



For TNOs, since they are assumed massless, the equation of 
motion is given by 



r n p 
h +&m s — = ^0mj 



v r,- 



■r k \ 



(3) 



where is the vector position of the TNO. By adopting the 
above equations we assumed that planets and TNOs are only per- 
turbed by the remaining planets, i.e., the TNOs are considered as 
test particles. The only exception will be the Pluto-Charon sys- 
tem barycenter, which will be considered as a planet because of 
its important mass. Indeed, some Plutinos may be pushed out of 
the 3:2 resonance by Pluto into close encou nters with Neptune 
dYu & Tremaindll999tlNesvornv et alj|2000h . 

In our simulations we will use the symplectic integrator from 
lLaskar & Robutell d2001t) . with an integration step-size of 0. 1 yr. 
The system is composed of 5 planets (Table lA"l2l ). 6 Trojans 
(TabHOt and 98 Plutinos (Table[Ql>. 

3.1. Resonant Motion 

Since both Trojans and Plutinos are resonant objects we will 
briefly recall here the b ases of resonant motion, following 
Murra y & Dermotj <|l999). Consider then a TNO in some res- 
onance with Neptune. We assume for simplicity, that Neptune is 



in a circular orbit and that all motion takes place in the plane of 
its orbit. We also ignore any perturbations between the two bod- 
ies as we are only interested in how resonant relationships lead 
to repeated encounters. 

We can examine the geometry of resonance for a general case 
by first considering two bodies moving around the Sun, in circu- 
lar and coplanar orbits. So, let us assume that 



n 
n 



p + q 



(4) 



where n and n ' are the mean motions of Neptune and the TNO, 
respectively (for Trojans p — 1 and q — 0, while for Plutinos 
p — 2 and q = 1). If the two bodies are in conjunction at time 
t — 0, the next conjunction will occur when (n - n')t = 2n, and 
the period, T co „, between successive conjunctions is given by 



In 



p t , = p + q T 



q 



(5) 



where T and T are the orbital periods of the two bodies. 

Now consider the case when e — 0, e' + 0, and m' + 0, 
where e denotes the eccentricity for Neptune and e ' and m' de- 
notes the eccentricity and longitude of pericentre of the TNO 
respectively. If the resonant relation 



{p + q)n' — pn — qm' — 

is satisfied, then we can rewrite Eq.© as 



n — m 



m P + q 



(6) 



(7) 



where n' - m' and n - m' are relative motions. These can be 
considered as the mean motions in a reference frame, co-rotating 
with the pericentre of the TNO. From the point of view of this 
reference frame, the orbit of the TNO is fixed or stationary. If 
the resonant relation given in Eq.© holds, the corresponding 
resonant argument is 



ip — (p + q)A' — pA — qm' 



(8) 



where A and A' denote the mean longitude of Neptune and the 
TNO, respectively. At a conjunction of the two bodies, A - A' 
and we have 



tp — q(A' — m') — q(A — ru') . 



(9) 



Thus, <p is a measure of the displacement of the longitude of con- 
junction from pericentre of the TNO. By computing the deriva- 
tive of the resonant angle tp, we get 



Cp = (p + q)n - pn — qm . 



(10) 



and ip — from Eq.©. In a more general situation, we will have 
<p + 0, but in order to preserve the resonant equilibrium cp will 
librate around an equilibrium position <po, obtained when ip — 0. 
The libration amplitude Atp will depend on the initial conditions 
and perturbations from the other bodies in the system and may 
reach large values. As a consequence, it is possible that the or- 
bits of two distinct bodies librating around different equilibrium 
positions intercept at some point. 

Along with the orbital parameters of the Trojans and Plutinos 
that we used in our simulations, in Table |A.3l and Table |A.4l re- 
spectively, we provide the equilibrium libration angle, the main 
libration period and amplitude of each TNO obtained over the 
next 250 Kyr. 
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3.2. Trojans 

Neptune Trojans are resonant objects in a 1:1 mean motion res- 
onance with Neptune. In our model we computed the motion of 
the six Neptune Trojans listed in Table lA.3l In Fig. (0, we show 
the behavior of all known Trojans, along time, in a co-rotating 
frame with Neptune for 100 Myr. Each dot shows the position of 
the TNO every 10 kyr. 

As expected, we see in Fig[3]that all Trojans orbit around 
the Lagrangian point L\, and execute tadpole-type orbits. This 
kind of orbits represents stable osc illations in the vicinity of th e 
Lagrangian equilibrium point (e.g. Giuliat ti Winter et alj |2007). 
The differences between the shape of their orbits, depend on the 
libration amplitude, but also on their orbital eccentricity and in- 
clination values. Trojan 2007VL305, that execute the most scat- 
tered orbit, also presents the largest eccentricity and inclination 
(<? = 0.062 and i = 28°). On the other hand, for small values 
of these two orbital parameters, the TNOs remain roughly in the 
path of Neptune's orbit, only changing its relative position to the 
planet due to libration. 



In Table lA.3l we provide the libration amplitude and period 
for all Trojans. While amplitudes can vary from only 6° to 26°, 
the periods of libration remain around 9 Kyr for all objects. 
Trojan 2001QR322 presents the largest libration amplitude and 
therefore moves further away of the equilibrium point L\. As a 
consequence, its orbit will be more susceptible of being desta- 
bilized by gravitational perturbations from the planets and other 
bodies in the system. Indeed, in our long-term numerical simu- 
lations (Sect. 14. Il l this TNO will abandon the Trojan orbit after 
1 12 Myr and become a Kuiper belt object. 

3.3. Plutinos 

Plutinos are resonant TNOs in a 3:2 mean motion resonance with 
Neptune. Thus, like Trojans, although they can cross the orbit of 
Neptune, they are protected from possible encounters with this 
planet. In Fig|4]we have drawn the typical path of a Plutino in 
the co-rotating frame of Neptune, for three different values of 
the eccentricity (e = 0.1, 0.2 and 0.3). 
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Fig. 3: Orbital evolution of the Neptune Trojans (in green) listed in 
Table lA.3l over 100 Myr in the co-rotating frame of Neptune (in blue). 
Each panel shows the projection of the TNO position every 10 Kyr in 
the orbital plane of Neptune, x and y are spatial coordinates centered 
in the star and rotating with Neptune, normalized by the Neptune-Sun 
distance. All Trojans orbit around the Lagrangian point L4, and execute 
tadpole-type orbits. More scattered orbits correspond to higher values 
of the eccentricity and inclination, while distance to the L 4 point depend 
on the libration amplitude. 



Fig. 4: Typical path of a Plutino (dotted line) in the rotating frame of 
Neptune (full line) for different eccentricity values (e = 0.1, 0.2 and 
0.3). The position of the Plutino was drawn for equal time intervals. 
Only high eccentricity values (e > 0.2) allow the Plutino to cross the 
orbit of Neptune. Due to the 3:2 mean motion resonance the trajectories 
are repeated every two orbits of the TNO around the Sun. 



These plots are drawn assuming that the Plutino is at exact 
resonance (ip = 0), which is not true, because the orbit is librat- 
ing around an equilibrium position <^o (Eq(8}. As a consequence, 
in a more realistic situation we will observe an oscillation of 
those paths as the one represented in Fig|5] In Table lA.4l we pro- 
vide the libration amplitude and period for all Plutinos. The equi- 
librium libration angle for them all is (po = +180°, but the am- 
plitudes of libration can be as small as 7° for Plutino 1996TQ66 
or as wide as 120° for Plutinos 1995QY9 and 2001KN77. The 
libration periods, Pm,, vary between 14.5 and 26.6 Kyr, the av- 
erage being around 20 Kyr. For comparison, the values for Pluto 
are A<p = 19.1° and P, ih = 19.9 Kyr. 



Plutino 




Fig. 5: Libration motion of the orbit of a Plutino. 
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In our model we computed the motion of about 100 Plutinos, 
whose orbital parameters are listed in Table |A.4l All objects 
present moderate eccent ricities and inclin ations, (e ~ 0.23 and 
i ~ 10.4°). According to Malhotraj d 1995b these values can be a 
consequence of the resonant mechanism of capture, during the 
residual planetesimal cleaning in the vicinity of the young giant 
planets. Due to Neptune's migration, the eccentricity and incli- 
nation of the Plutinos are pumped after capture in resonance. 

In Fig|6]we show the behavior of six Plutinos, along time, in 
a co-rotating frame with Neptune for 100 Myr, each dot show- 
ing the position of the object every lOKyr. Because Plutinos 
are much more numerous than Trojans, we can only represent a 
small fraction of them. We have chosen the most representative 
cases, namely, the Plutino with the smallest and widest libra- 
tion (1996TP66 and 2001KN77, respectively), the Plutinos with 
smallest and highest eccentricity (2003VS2 and 2005GE187, re- 
spectively), and the Plutinos with smallest and highest inclina- 
tion (2002VX130 and 2005TV189, respectively). 
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Fig. 6: Orbital evolution of some Plutinos (in red) taken from Table lA.4l 
over 100 Myr in the co-rotating frame of Neptune (in blue). Each panel 
shows the projection of the TNO position every 10 Kyr in the or- 
bital plane of Neptune, x and y are spatial coordinates centered on the 
star and rotating with Neptune, normalized by the Neptune-Sun dis- 
tance. We have chosen the most representative cases, namely, Plutino 
1996TP66, with the smallest libration amplitude (Aip = 1.2°, e = 0.328, 
i = 5.6°), Plutino 2001KN77, with the widest libration amplitude 
(Atp = 120.4°, e = 0.242, i = 2.4°), Plutino 2003VS2, with the small- 
est eccentricity (e = 0.072, i = 14.8°), Plutino 2005GE187, with the 
highest eccentricity (e = 0.329, i = 18.2°), Plutino 2002VX130, with 
the smallest inclination (e = 0.220, i = 1.3°), and Plutino 2005TV189, 
with the highest inclination (e = 0.186, ( = 34.5°). 



As expected, depending on the eccentricity values, the orbits 
of the Plutinos are all in good agreement with the paths shown 
in Fig0] Because of the libration motion of the orbits (FigO 
their trajectories approach or cross the Lagrangian points L4 and 
L5, that is, the Plutinos orbits share the same spatial zone as the 
Neptune Trojans. 



4. Numerical Simulations 

In Sectf2] we have shown that, although Neptune Trojans and 
Plutinos appear to have different origins, some of their proper- 
ties indicate that the two types are quite identical, suggesting 
that there must be some sort of communication between them. 
Indeed, in the previous section we saw that due to the libra- 
tion motion of the orbits there is a wide zone of spatial overlap 
around the Lagrangian point L4 of Neptune. As a consequence, 
we may expect close encounters and collisions between the two 
kinds of TNOs to occur at a higher rate than in the remaining 
Kuiper belt, resulting in a change of the size distributions of the 
two populations. 

In order to test this possibility we simulated the long-term 
future evolution of the outer Solar System for 1 Gyr. The orbits 
of the outer planets, the Neptune Trojans and the Plutinos are 
integrated simultaneously according to the model described in 
the beginning of Sectj3] 

4.1. Stability of the Neptune Trojans 

The stability of the Neptune Trojans orbits is an important is- 
sue on the dynamic s of the outer Solar System. According to 
iDvorak et al.l d2007l) Trojans with low-inclined orbits are less 
stable. The stability area around L4 and L5 disappears after 
about 10 8 yr for low inclinations, while this stability zone is still 
present for about 10 9 yr for large inclinations. More precisely, 
it was concluded that there exists a region (20° < i < 50°) of 
higher stability for the Neptune Trojans, although only two have 
presently been found in this region (Table lTOT i. 

During our numerical simulations all Trojans remained sta- 
ble within the limits shown in Figf3]during 1 Gyr, except Trojan 
2001QR322, which escapes from the Lagrangian L4 point af- 
ter about 112 Myr. This las t obse rvation was somehow un- 
expected, since Chiang et al.l d2003l) concluded that this same 
Neptune Trojan was stable over 1 Gyr. This difference of behav- 
iors is a consequence of a modification in the initial conditions 
(Table fATJt , obtained with new observational data. This event is 
more or less in agreement with the results from IDvorak et al.l 
(12.007b . because Trojan 2001QR322 is the one with smallest in- 
clination (z = 1.3°). Although other Trojans with identical incli- 
nation values remained stable, we notice that Trojan 2001QR322 
also has the widest libration amplitude {A<p = 26°) and its or- 
bit is therefore more susceptible of being disturbed by planetary 
perturbations. This suggests that the stability of the Trojans is 
smaller for low-inclined orbits, but also for large libration am- 
plitudes. 

In FigQwe plot the long-term evolution of the orbital pe- 
riod and the eccentricity of the Trojan 2001QR322. The behav- 
ior of this object is extremely regular, until it undergoes a sud- 
den increment of the eccentricity, which removes it form the 1 : 1 
mean motion resonance with Neptune. Interestingly, just after 
escaping the resonance, the semi-major axis of this TNO is tem- 
porarily stabilized near the 3:2 resonance with Neptune, i.e., its 
orbit becomes very close to the Plutino's orbits. However, con- 
trary to regular Plutinos, the eccentricity undergoes important 
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chaotic variations from nearly zero up to 0.3. This regime lasts 
for about 300 Myr, time after which the eccentricity grows to 
almost 0.8. As a consequence, close encounters with the plan- 
ets become possible and the TNO leaves the area near to the 
3:2 mean motion resonance. Later on, the same TNO seems to 
be captured in a 9/2 mean motion resonance with Neptune, and 
stays there for about 100 Myr. Finally, the orbit is again destabi- 
lized and the eccentricity grows to extreme values, and the TNO 
turns into a comet. 



Table 1 : List of all unstable Plutinos during 1 Gyr. 



Plutino 



1e+08 2e+08 3e+0S 

time/yr 







20O1QR322 , 






! 






i 









1e+08 2e+08 3e+08 4e+08 5e+08 

time/yr 

Fig. 7: Long-term evolution of the orbital period (over Neptune's or- 
bital period) and the eccentricity of Trojan 2001QR322 for 525 Myr. 
Initially a Trojan, the orbit of this TNO is not stable and quits the 1:1 
mean motion resonance after about 112 Myr. Then it jumps to different 
configurations until it becomes a comet and eventually collides with a 
planet or the Sun. In this simulation it is ejected from the Solar System 
after 525 Myr. 



4.2. Stability of the Plutinos 

The orbits of the Plutinos are believed to be stable as long as 
the lib ration amplitudes are smaller than 120° dNesvornv & Roia 
2000). For larger libration widths, the planetary perturbations 
will remove the TNO from its orbit in short time intervals. The 
same phenomena was observed for Trojans. 

During our numerical simulations over 1 Gyr, 10 Plutinos 
over 98 quit their orbits (about 10%). In Table[T]we list all these 
bodies, as well as their libration width (Table lA~4l ). Contrary to 
expectations, we observe that only two Plutinos present libration 



time (Myr) 



a<pO 



2000FV53 


69 


0.168 


117.76 


2004FU148 


230 


0.235 


94.67 


2002GE32 


500 


0.232 


73.14 


1998WZ31 


590 


0.165 


76.61 


2004EW95 


595 


0.320 


53.31 


1995QY9 


680 


0.262 


120.22 


2003TH58 


710 


0.088 


65.29 


1993SB 


720 


0.317 


53.19 


2002XV93 


820 


0.127 


42.70 


2003UT292 


840 


0.292 


80.69 



amplitudes around 120°. The minimum libration width observed 
is around 40°. A possible explanation is t hat because Plutinos are 
not alone in their orbits an d, according to lYu & Tremaind (1999) 
and iNesvornv et al.l d2000l) . some Plutinos may be pushed out 
of the resonance by Plu to into close encounters with Neptune. 
Indeed, stability maps bv lNesvornv & R oig (2000) only take into 
account the effect of the four giant planets. 

In FigJS] we plot the long-term evolution of the orbital pe- 
riod and the eccentricity of TNO 2000FV53. The initial 3:2 res- 
onant configuration is abandoned just after 69 Myr, time after 
which the eccentricity increases progressively. The TNO then 
undergoes close encounters with the planets, which will increase 
the eccentricity even more, until it can reaches very high val- 
ues. The Plutino 2000FV53 then becomes a short-period comet 
and eventually collides with a planet or the Sun. This mecha- 
nism has already been described as the responsible for the pro- 
vision of short period comets into the inner Solar System (e.g. 
Morbidelli 1997). In our simulation the aforementioned TNO is 
ejected from the Solar System after 150 Myr. 

4.3. Orbital overlap between Trojans and Plutinos 



In Sect l3.3l we saw that because of libration, the Plutinos' orbits 
can approach the Lagrangian point L\. In order to check the ex- 
tent of the orbital overlap between Neptune Trojans and Plutinos, 
in Fig|9]we plotted simultaneously their orbits in a co-rotating 
frame with Neptune. We used the Trojan 2007VL305 for all rep- 
resentations, since it has the most scattered orbit, maximizing the 
possibility of orbital merging with a Plutino. For the Plutinos we 
used the same bodies as in Sect l3.3l which correspond to the ex- 
treme cases of libration, eccentricity, and inclination. The only 
exception is that Plutino 1996TP66 (with the smallest libration 
width) has been replaced by Pluto, since a small libration ampli- 
tude does not allow the Plutino to flyby the Lagrangian points. 

Since the TNOs are not necessarily in the same orbital plane 
as Neptune (especially for those having large inclination val- 
ues), two types of plots have been made: one where we plot- 
ted the projection of the TNO position in the orbital plane of 
Neptune (Fig(9^), and another where we plotted the projection 
of the TNO in the orbital plane of Neptune, only when the dis- 
tance to this plane is smaller then 10~ 3 AU (Fig|9}3), that is, less 
than 150 000 km, half of the Earth-Moon distance. Indeed, most 
of our Trojans have low inclinations and therefore lie close to 
Neptune's orbital plane most of the time. As a consequence, near 
this plane close encounters with Plutinos are maximized. 

The importance of plotting the two situations is clearly illus- 
trated by the behavior of Pluto (FigO. At first glance, looking 
only to the projection on Neptune's orbital plane, we observe 
a large zone shared by the orbits of the two kinds of TNOs, 
suggesting that close encounters may be a regular possibility. 
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Fig. 9: Orbital evolution of the Trojan 2007VL305 (in green) and Pluto plus most of the Plutinos in Fig[6](in red), over 1 Gyr in the co-rotating 
frame of Neptune (in blue). Left panel (a) shows the projection of the TNO position every 100 Kyr in the orbital plane of Neptune, while right 
panels (b) shows the projection of the TNO in the orbital plane of Neptune every 10 Kyr, only when the distance to this plane is smaller then 
10 -3 AU. We observe that orbital overlap is favored for Plutinos with high libration amplitudes, high eccentricity values and low inclined orbits. 



However, when we restrain the plot only to Neptune's orbital 
plane, we observe that Pluto is never on this plane when it ap- 
proaches L4, preventing any close encounter with Trojans. 

From the analysis of Figj9j we can also conclude that 
Plutinos with large libration amplitudes (e.g. 2001KN77) maxi- 
mize the chances of intercepting Trojans. This behavior was ex- 
pected, as their orbits invade a large zone around the Lagrangian 
point L\. 

On the other hand, Plutinos with low eccentricity values (e ~ 
0.1; e.g. 2003VS2), always avoid Trojans independently of its 
libration amplitude or orbital inclination, since the small eccen- 
tricity prevents them from crossing Neptune's orbit. According 
to FiglU this result could also be expected, because for low ec- 
centricity there is no interception of the Plutino and Neptune's 
orbits. 

Finally, when comparing the behavior of Plutinos in low- 
inclined orbits (2001KN77 and 2002VX130) with high-inclined 
ones (Pluto, 2003 VS2, 2005GE187 and 2005TV189) we con- 
clude that large orbital inclination decreases the chances of close 
encounters because the Plutino is never close to Neptune's or- 
bital plane when it crosses the orbit of the planet. The fact that 
Plutinos with low inclination share the Trojan space was ex- 
pectable, as both TNOs remain close to the same orbital plane. 
However, Plutinos with high inclination could also approach the 
orbital plane of Neptune at the moment they are close to the 
Lagrangian points, which is not really observed. 



From the above analysis (and also for the remaining 5 
Trojans and 92 Plutinos), we empirically conclude that orbital 
overlap between Trojans and Plutinos is favored for Plutinos 
with high libration amplitudes, high eccentricity values and low 
inclined orbits. 

4.4. Collisions between Neptune Trojans and Plutinos 

In order to directly check if close encounters between the 
Neptune Trojans and Plutinos can occur, and how often during 
1 Gyr of numerical simulations, we computed the distance be- 
tween all bodies after each step-size. 

For that purpose we arbitrarily selected two critical dis- 
tances, one d\ < 2 x 10~ 5 AU (~ 3 000 km), for which we 
assume that the two TNOs effectively collide, and a second 
d 2 < 2 x 10 3 AU (~ 300000 km), for which the two bodies 
do not collide, but become closer than the Earth-Moon distance. 
Assuming a body with the mass density of water ice located at 
Neptune's distance to the Sun, d\ and &i correspond approxi- 
mately to the Hill radius of an object with diameter of 1 .5 km 
and 150 km, respectively. Alternatively, d 2 is equivalent to the 
Hill radius of a body with 0.01 % of Pluto's mass at the same dis- 
tance from the Sun as Neptune. TNOs with diameters of 150 km 
or larger when approaching at distances less than d 2 will be sig- 
nificantly perturbed by their mutual gravity and our model de- 
scribed in Sectj3] will no longer apply. We assume that TNOs 
undergoing such close encounters may effectively collide, or de- 
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Table 2: Close encounters between TNOs during 1 Gyr (d mj „ < d 2 ). 



Type Body 1 


Body 2 


time 








(Gyr) 


(km) 



2e+07 4e+07 6e+07 8e+07 

time/yr 



1e+08 1.2e+08 1.4e+08 




2e+07 4e+07 6e+07 8e+07 1e+08 1.2e+08 1.46+08 
time/yr 

Fig. 8: Long-term evolution of the orbital period (over Neptune's orbital 
period) and the eccentricity for Plutino 2000FV53 for 150 Myr. Initially 
a Plutino, the orbit of this TNO is not stable and quits the 3:2 mean mo- 
tion resonance just after about 69 Myr. Its eccentricity then increases 
progressively and the TNO becomes a short-period comet that eventu- 
ally collides with a planet or the Sun. In this simulation it is ejected 
from the Solar System after 150 Myr. 



viate considerably from their initial orbits and quit the resonant 
configuration with Neptune. 

After 1 Gyr of simulations we did not observe any event 
for which the minimal distance between two bodies, d min , is 
lower than d\, and registered only 14 close encounters for which 
d,„i„ < di- The results are listed in Table|2] However, these re- 
sults cannot be seen as definitive, but rather as minimal estima- 
tions of close encounters. Indeed, since our step-size is 0.1 yr, in 
a circular orbit a Trojan will travel about 0.1 AU per step-size. 
As a consequence, two TNOs may effectively collide between 
two step-sizes and our program is unable to detect it. The results 
listed in Table|2]must then be seen as indicative of the possibility 
of collisions and not as conclusive. 

Assuming a constant speed for TNOs and a uniform distri- 
bution of their relative minimal distances, we roughly estimate 
the real number of close encounters d mm < di to be 50 times 
more frequent than those listed in Table|2] (0.1 AU/c/2 = 50). 
Effective collisions (af m ,„ < d\) should also be more frequent in 
the same proportion. Thus, since d\ = \0~ 2 d2, the results showed 
in Table|2]for d mm < di can be seen as a rough indicator of the 
real number of effective collisions occurring between Neptune 
Trojans and Plutinos. 



T-P 


2005TN53 


2002VU130 


0.244 


270 226 


P-P 


2003SO317 


1996SZ4 


0.278 


142 977 


T-P 


2005TO74 


1995HM5 


0.298 
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P-P 


1998HH151 


2001KX76 


0.320 


196 351 


P-P 


2001QH298 


1993SC 


0.406 


250 561 


P-P 


1998HH151 


2003AZ84 


0.455 


232 759 


P-P 


2001KB77 
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0.494 


281 273 


P-P 


2003SR317 


1994TB 


0.519 


271635 


P-P 


2001KQ77 


2005TV189 


0.528 


193 516 


P-P 


1998WS31 


2004FW164 


0.542 
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P-P 


1998WV31 


1998HK151 


0.588 


236512 


T-T 


2004UP10 


2005TO74 


0.849 


272 851 


T-T 


2004UP10 


2006RJ103 


0.864 


159 886 


P-P 


2004EH96 


1993SC 


0.986 


226 779 



Among the 14 "collisions" listed in Table[2] two were be- 
tween Trojans, two between a Trojan and a Plutino, and the re- 
maining ten between Plutinos. It is not a surprise that Trojans or 
Plutinos also undergo close encounters between each other, be- 
cause of the libration of their orbits (Pig.|5j. The two Plutinos 
that encounter Trojans are 1995HM5 (e = 0.26, i = 4.8°, 
Atp = 70°) and 2002VU130 (e = 0.21, i = 1.4°, Aip = 116°), 
both having small orbital inclinations and large values for the 
eccentricity and libration amplitude. This confirms our predic- 
tions from Sect l4.3l that collisions between Trojans and Plutinos 
are favored for Plutinos with high libration amplitudes, high ec- 
centricity values and low inclined orbits. 

The fact that we count more close encounters between two 
Plutinos than between a Trojan and a Plutino, cannot be seen 
as an indicator that this last kind of encounter is less frequent. 
Indeed, in our simulations the number of Trojans (6) is about 16 
times smaller than the number of Plutinos (98). Therefore, there 
are roughly 16 times more chances of observing an encounter be- 
tween two Plutinos. If our simulations had as much Trojans as it 
has Plutinos, we could then probably expect to observe about 32 
Trojan-Plutino encounters. As a consequence, from the results 
listed in Table|2]we infer that this kind of encounter is roughly 
3 times more frequent than a Plutino-Plutino one. The same ap- 
plies to the encounters between Trojans. In our simulation they 
are also about 16 times less probable than a Trojan-Plutino en- 
counter and we infer that they should be about 50 times more fre- 
quent than encounters between Plutinos. Only by using a model 
with an identical initial number of Trojan and Plutinos, would 
allow us to determine the exact proportions for each kind of en- 
counter. 



5. Conclusions 

In this work we aimed to verify if both Plutino and Neptune 
Trojan populations can be merged, and how frequently does that 
occur, and examine if there could be a connection between such 
possibility and the observed properties of the two (sub-) popu- 
lations of TNOs in question. We analyzed the available colors 
and absolute magnitudes of Plutinos and Neptune Trojans. The 
nonexistence of significant albedo diversity among the large ma- 
jority of these objects is assumed, hence we interpreted absolute 
magnitudes as an equivalent to to object size. We find that: 

(i) there are no intrinsically bright (large) Plutinos at small in- 
clinations; 
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(ii) there is an apparent excess of blue and intrinsically faint 
(small) Plutinos; 

(iii) Neptune Trojans possess the same blue colors as Plutinos 
within the same (estimated) size range do. 

From these results, we have hypothesized that there might 
have been some strong collisional interaction between Neptune 
Trojans and (1) the presently small Plutinos — with the assump- 
tion that such collisions created only blue colored objects — , or 
(2) the low inclined Plutinos — with the assumption that such 
collisions created both blue and red objects. Paradoxically, the 
previous two opposite assumptions both seem to have observa- 
tion support among other TNOs. 

In order to differentiate between these two scenarios we per- 
formed a numerical simulation over 1 Gyr of the future evolution 
of the outer Solar System composed by 5 planets, 6 Neptune 
Trojans, and 98 Plutinos. 

By plotting simultaneously the orbits of the Neptune Trojans 
and Plutinos in a co-rotating frame with Neptune, it becomes 
clear that there is a large overlap of the orbits of the two kinds 
of TNOs. A more detailed analysis revealed, however, that close 
encounters with the Neptune Trojans are favored for Plutinos 
with high libration amplitudes, high eccentricity values, and low 
inclined orbits. Tables |A~3l and |A.4l show the libration amplitudes 
and periods computed for these objects. 

After 1 Gyr of numerical simulations we registered 2 Trojan- 
Plutino, 2 Trojan-Trojan, and 10 Plutino-Plutino close ap- 
proaches, i.e. less than ~ 300000 km. Since we have used the 
currently known objects, Trojans are much less numerous than 
Plutinos in our simulation. If the number of objects among each 
population is similar then we roughly infer that Trojan-Plutino 
and Trojan-Trojan "collisions" can be about three and fifty times 
more frequent, respectively, than Plutino-Plutino collisions. 

The collision rates between Neptune Trojans and between 
Neptune Trojans and Plutinos might explain both why we ob- 
serve a small number of Neptune Trojans and why there is an 
absence of large Neptune Trojans: a strong collisional evolution 
possibly played an important role by shattering and/or depleting 
Neptune Trojans. 

Our results also show that Plutinos in low-inclined orbits 
have more chances of colliding with Neptune Trojans. This 
result gives no support for a (mutual) collisional origin of 
the equal-sized and equal-colored Neptune Trojans and small 
Plutinos, as the latter are equally spread in inclination (see Sect. 
|2j scenario #1). On the other hand, it gives plausibility to the ori- 
gin of the concentration of small Plutinos with i < 8 - 13° as be- 
ing a consequence of some collisional interaction with Neptune 
Trojans(see Sect.ffJ scenario #2). 

During our numerical simulations we also observed that the 
orbit of Trojan 2001QR322 became unstable as well as the orbits 
of 10 more Plutinos. The stability of Ne ptune Trojans appe ars to 
be enhanced for high inclination values dDvorak et al.l2007b . and 
also for small libration amplitudes. On the other hand, Plutinos 
seem to be essentially pushed out of their resonan c e by P luto, 
in conformity with th e results of lYu & Tremaind d 19991) and 
iNesvornv et ai](l2000l) . 

This work's results were derived for objects considered as 
test particles (except for Pluto) and the number of "collisions" 
were extrapolated from the amount of close encounters after 
each integration step-size. Accurate estimations for the amount 
of collisions between Neptune Trojans and Plutinos can only be 
made with the inclusion of mutual gravitational interactions be- 
tween all objects, and using the real spatial and size distributions 
of these objects. 



Nonetheless, our sketchy analysis indicates that under cer- 
tain assumptions, which have parallel in what has been inferred 
for other TNOs, shattering collisions involving Neptune Trojans 
might have played a crucial role on the creation of the size- 
inclination asymmetries observed among Plutinos. We recall that 
we have disregarded the possibility that these asymmetries might 
have been created by collisional interaction with some other 
sub-population of TNO s than the Neptune Trojans in vie w of 
the results obtai ned bv lThebault & Dor essoundiram (120031) and 
Thebault (2003). If those works come to be revealed as inade- 
quate approximations of the relative collision rates suffered by 
TNOs our inference on a possible cause for the size-inclination 
distribution of Plutinos looses its ground. Further, if Neptune 
Trojans and low inclined Plutinos do not possess identical spin 
rate versus size distributions, which should be distinct from the 
high er inclined Plutinos, t hen our suggestions cannot hold either 
(e.g. lFarinella et al"1ll98ll) . More detailed studies on the interac- 
tion between Neptune Trojans and Plutinos should be attempted. 
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Table A. 2: Orbital data for the Giant Planets and Pluto at JD 2454200.50 (http://ssd.jpl.nasa.gov/horizons.cgi I. 



Name 


a(AU) 


e 


i (deg) 


M (deg) 


u (deg) 


a (deg) 


m (M Q ) 


Jupiter 


5.20219308 


0.04891224 


1.30376425 


240.35086842 


274.15634048 


100.50994468 


9.5479194xl0~ 4 


Saturn 


9.54531447 


0.05409072 


2.48750693 


45.76754755 


339.60245769 


113.63306105 


2.8586434xl0~ 4 


Uranus 


19.19247127 


0.04723911 


0.77193683 


171.41809349 


98.79773610 


73.98592654 


4.3558485xl0- 3 


Neptune 


30.13430686 


0.00734566 


1.77045595 


293.26102612 


255.50375800 


131.78208581 


5.1681860x10^ 


Pluto 


39.80661969 


0.25440229 


17.121129 


24.680638 


114.393972 


110.324800 


6.5607561xl0- lJ 



Table A.3: Orbital data for the Neptune Trojans at JD 2454200.50 (ftp://ftp.lowell.edu/pub/elgb/astorb.html I. 



# 


Name 


a(AU) 


e 


i (deg) 


M (deg) 


oj (deg) 


Q. (deg) 


Pm (Kyr) 


<Po (deg) 


A<p (deg) 


1 


2001QR322 


30.190 


0.029 


1.3 


60.2 


154.8 


151.7 


9.23 


68.10 


25.90 


2 


2004UP10 


30.099 


0.025 


1.4 


334.1 


2.2 


34.8 


8.86 


61 .44 


10.5 


3 


2005TN53 


30.070 


0.062 


25.0 


280.3 


88.6 


9.3 


9.42 


58.95 


6.61 


4 


2005TO74 


30.078 


0.051 


5.3 


260.1 


306.9 


169.4 


8.80 


60.91 


6.88 


5 


2006RJ103 


29.973 


0.028 


8.2 


226.6 


35.4 


120.8 


8.87 


60.45 


6.13 


6 


2007VL305 


29.956 


0.061 


28.1 


348.5 


216.1 


188.6 


9.57 


61.08 


14.26 



AJ.C. Almeida et al.: Neptune Trojans and Plutinos: colors, sizes, dynamics, and collisions 
Table A.4. Orbital data for the Plutinos at JD 2454200.50 (ftp://ftp.lowell.edu/pub/elgb/astorb.html I. 
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52 


2001VN71 


39.287 


0.243 


18.692 


359.265 


1.438 


70.441 


22.76 


179.14 


49.91 


53 


2001YJ140 


39.282 


0.290 


5.980 


358.246 


129.452 


319.434 


19.73 


179.34 


74.73 


54 


2002CE251 


39.543 


0.272 


9.294 


347.383 


215.444 


342.554 


14.46 


-177.96 


103.83 


55 


2002CW224 


39.185 


0.243 


5.668 


293.038 


156.088 


1.759 


21.50 


180.41 


43.34 


56 


2002GE32 


39.569 


0.232 


15.670 


287.290 


103.320 


203.739 


20.16 


-180.41 


73.14 


57 


2002GF32 


39.497 


0.172 


2.779 


111.703 


54.825 


44.317 


18.89 


-181.23 


91.71 
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Table A.4 - Continuation from previous page 



# 


Name 


a (AU) 


e 


i (deg) 


M (deg) 


to (deg) 


O (deg) 


Pub (Kyr) 


<Po (deg) 


Aip (deg) 


58 


2002GL32 


39.723 


0.131 


7.070 


4.155 


192.923 


11.080 


21.80 


-179.71 


53.86 


59 


2002GV32 


39.797 


0.198 


5.373 


349.937 


173.000 


79.151 


21.32 


-178.84 


43.97 


60 


2002GW31 


39.413 


0.239 


2.640 


87.855 


198.315 


227.267 


19.50 


179.04 


80.27 


61 


2002GY32 


39.716 


0.095 


1.799 


16.554 


337.273 


225.560 


22.49 


-179.65 


24.46 


62 


2002VD138 


39.403 


0.151 


2.784 


41.170 


36.198 


315.079 


20.18 


178.95 


76.03 


63 


2002VE95 


39.132 


0.285 


16.346 


8.236 


206.775 


199.855 


21.38 


181.71 


56.82 


64 


2002VR128 


39.313 


0.265 


14.035 


60.630 


287.630 


23.108 


22.00 


177.26 


18.06 


65 


2002VU130 


39.022 


0.211 


1.373 


258.236 


281.602 


267.864 


16.27 


181.94 


115.85 


66 


2002VX130 


39.325 


0.220 


1.322 


359.559 


106.036 


296.778 


21.18 


178.52 


46.47 


67 


2002XV93 


39.204 


0.127 


13.286 


267.105 


165.694 


19.121 


21.96 


176.46 


42.70 


68 


2003 AZ84 


39.413 


0.181 


13.596 


215.203 


15.040 


252.143 


22.84 


179.64 


44.63 


69 


2003FB128 


39.821 


0.260 


8.867 


38.023 


306.210 


209.482 


18.96 


-179.28 


88.10 


70 


2003FF128 


39.831 


0.221 


1.911 


333.787 


169.763 


91.731 


20.25 


-179.16 


65.21 


71 


2003FL127 


39.337 


0.233 


3.500 


146.187 


55.981 


314.317 


20.42 


178.97 


63.50 


72 


2003HA57 


39.648 


0.176 


27.584 


0.276 


5.641 


199.708 


25.80 


-186.74 


48.17 


73 


2003HD57 


39.697 


0.183 


5.612 


22.980 


137.734 


34.397 


21.15 


-179.59 


53.14 


74 


2003HF57 


39.619 


0.198 


1.422 


23.776 


123.869 


48.151 


21.10 


-178.70 


50.19 


75 


2003QB91 


39.210 


0.197 


6.493 


132.226 


80.949 


136.756 


19.18 


182.13 


85.39 


76 


2003QH91 


39.540 


0.152 


3.652 


108.349 


267.910 


286.690 


18.07 


182.40 


105.60 


77 


2003QX111 


39.403 


0.146 


9.534 


85.767 


101.157 


157.498 


21.42 


183.43 


44.50 


78 


2003SO317 


39.318 


0.276 


6.573 


42.879 


111.721 


187.304 


17.98 


179.61 


93.04 


79 


2003SR317 


39.426 


0.168 


8.357 


50.834 


117.872 


175.647 


19.69 


181.16 


79.06 


80 


2003TH58 


39.240 


0.088 


27.994 


15.262 


166.733 


251.400 


22.17 


175.75 


65.29 


81 


2003UT292 


39.095 


0.292 


17.570 


332.561 


256.085 


211.037 


19.57 


181.5 


80.69 


82 


2003UV292 


39.197 


0.214 


10.998 


46.554 


121.092 


235.689 


21.33 


177.93 


46.02 


83 


2003VS2 


39.266 


0.072 


14.800 


4.586 


112.370 


302.667 


25.86 


179.66 


27.96 


84 


2003WA191 


39.157 


0.232 


4.521 


13.619 


226.152 


179.282 


21.03 


179.04 


51.52 


85 


2003WU172 


39.058 


0.254 


4.148 


338.412 


101.101 


10.422 


17.72 


179.89 


104.89 


86 


2004DW 


39.300 


0.222 


20.594 


162.480 


72.412 


268.724 


21.07 


182.84 


67.54 


87 


2004EH96 


39.639 


0.283 


3.128 


11.376 


301.717 


226.249 


18.35 


-178.92 


89.89 


88 


2004EJ96 


39.795 


0.244 


9.327 


323.832 


231.714 


18.733 


18.53 


-178.74 


93.32 


89 


2004EW95 


39.670 


0.320 


29.241 


344.317 


204.373 


25.751 


22.40 


-176.50 


53.31 


90 


2004FU148 


39.867 


0.235 


16.623 


318.370 


81.079 


197.868 


18.82 


-178.18 


94.67 


91 


2004FW164 


39.752 


0.163 


9.099 


359.209 


7.466 


197.971 


20.74 


-178.92 


71.70 


92 


2005EZ296 


39.647 


0.155 


1.773 


322.078 


211.628 


24.433 


18.12 


-177.57 


103.09 


93 


2005EZ300 


39.721 


0.243 


10.319 


311.819 


252.031 


357.124 


16.54 


-179.18 


114.11 


94 


2005GA187 


39.653 


0.221 


18.714 


290.525 


281.974 


27.716 


21.60 


-182.86 


46.21 


95 


2005GB 187 


39.743 


0.242 


14.659 


6.049 


349.134 


217.084 


22.46 


-180.87 


35.64 


96 


2005GE187 


39.660 


0.329 


18.222 


324.766 


85.071 


205.378 


20.93 


-179.44 


38.86 


97 


2005GF187 


39.803 


0.262 


3.906 


337.514 


134.526 


128.952 


21.17 


-179.76 


38.06 


98 


2005TV189 


39.286 


0.186 


34.455 


359.741 


186.066 


246.224 


28.54 


164.83 


33.49 
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